Intrinsic vs. laboratory frame description of the deformed nucleus 48 Cr. 
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The collective yrast band of the nucleus 48 Cr is studied 
using the spherical shell model and the HFB method. Both 
approaches produce basically the same axially symmetric in- 
trinsic state up to the - accurately reproduced - observed 
backbending. Agreement between both calculations extends 
to most observables. The only significant discrepancy comes 
from the static moments of inertia and can be attributed to 
the need of a more refined treatment of pairing correlations 
in the HFB calculation. 
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The study of the collective behavior of deformed nuclei 
is a classical problem in Nuclear Physics. Traditionally, 
mean field descriptions in the intrinsic frame have been 
favoured, as they take naturally advantage of the spon- 
taneous breakdown of rotational symmetry. The price to 
pay for the gain in physical insight is the loss of angular 
momentum as good quantum number. 

In the laboratory frame description, as provided by 
spherical shell model calculations (SM), angular momen- 
tum is conserved but the physical insight, associated to 
the existence of an intrinsic state is lost, except in the 
very rare cases where Elliott's SU3 symmetry [Q oper- 
ates. Furthermore, the approach suffers form numeri- 
cal limitations. Hence, so far, it had been implemented 
mostly in regions such as the p and sd shells where the 
number of active particles is too small for collective fea- 
tures to become dominant. Nonetheless, there are a few 
nuclei - such as 20 Ne and 24 Mg - that are well repro- 
duced by the SM calculations and do exhibit collective 
properties, whose origin can be traced to the approximate 
validity of the SU3 symmetry, for which the relationship 
between the intrinsic and laboratory frame descriptions 
is well understood. 

In regions where the SU3 symmetry is poorly re- 
spected, as in in the pf shell g], the study of poten- 
tially good "rotors" was impaired by lack of experi- 
mental evidence, and by the difficulty of an exact SM 
treatment beyond 5 active particles. The situation has 
changed through recent measurements || demonstrating 
that 48 Cr is a good rotor up to spin J = 10 where the 
yrast band bends back. This behavior is reminiscent of 
the situation in much heavier deformed nuclei. Simul- 
taneously, full pf calculations B have become available, 
that reproduce in detail the observed properties of A=48 
isobars, and in particular those of 48 Cr. 

Therefore, this nucleus provides an unique testing 



ground to compare the SM (laboratory frame) descrip- 
tion of permanent deformation with Cranked Hartree- 
Fock- Bogoliubov (CHFB) calculations [5 with the finite 
range density dependent Gogny force [6|; which repre- 
sent the (self-consistent) state of the art formulation of 
the intrinsic frame approach. 

From the comparison it should be possible to obtain a 
better understanding of the intrinsic structure of the SM 
solutions, which in turn, may indicate in what sense the 
CHFB description falls short of an exact one. 

Computational procedures. In the Spherical Shell 
Model (SM) Cr is described in a 0hu> space, i.e. eight 
particles are allowed to occupy all the states available 
in the pf shell (1963461 states). The effective interac- 
tion is given by a minimally modified version of the Kuo- 
Brown's G-matrix H denoted KB3 in Q . The single par- 
ticle energies are taken from the Ca experimental spec- 
trum. The effect of core polarization on the quadrupole 
properties is taken into account by the use of effective 
charges = 1.5, q v = 0.5. The Hamiltonian is treated 
by the Lanczos method and diagonalized by the code 
ANTOINE g. 

In the intrinsic frame calculations we have used the Self 
Consistent Cranking Hartree- Fock- Bogoliubov method 
( CHFB) with the density dependent Gogny force. The 
CHFB equations determining the mean field intrinsic 
state \4>ui) are obtained by imposing the condition that 
the mean value of the Routhian be stationary against 
small variations of the intrinsic state, i.e., 



S{<t> u \H-wj x -\ N N-X z Z\<f> u ) = 0. 



(1) 



The Lagrange multipliers lj, A at and Xz are determined 
by the usual angular momentum and particle number 
constraints {^\J X \^) = + 1), {^\N\^) = N 

and = Z. 

The HFB wave functions have been expanded in a tri- 
axial harmonic oscillator basis I ) with different 

oscillator lengths. Ten oscillator shells are included in 
order to ensure the convergence of the mean field re- 
sults. The parameters of the Gogny force used in this 
calculation were adjusted more than ten years ago to re- 
produce ground state bulk properties of nuclei (DS1 set 
H). Without further changes, this force has proven ca- 
pable of describing successfully many phenomena, and in 
particular high spin behaviour j^j. 

In order to understand more qualitatively the physics 
involved and to make contact with the Shell Model cal- 
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culations we have computed the following quantities in a 
spherical representation of the basis: 
The "fractional shell occupancy" 

m=j 

= 27+T zZ (^Kljm^ljmlM, (2) 
m=-j 

the "shell contribution to (Jx)" 
3x(n,l,j) = 

^ (^)(nijm),(ni J W)(^|c+y m C ri ij m ^ w ) (3) 
m,m' 

and the "shell contribution to the quadrupole moment" 
Q2o(n,l,j;ri,l',j') = 

^ ( l ?2o)(niim),(n'i'i'm')^ w l C nijm C "' i 'j' m 'l^)- ^ 
m,m' 

In the above formulae \<j>J) is the intrinsic CHFB wave 
function expressed in the triaxial basis and c+y m are the 
operators creating a particle in the harmonic oscillator 
orbit \nljm) with oscillator length bo = (bxbybz) 1 ^ 3 . In 
order to obtain these quantities the triaxial basis has been 
expanded in a spherical one following a procedure similar 
to that of Ref. [lofl . As the triaxial basis has, in gen- 
eral, different oscillator lengths the expansion contains 
in principle an infinite number of terms. In our case, an 
efficient truncation is obtained by allowing the spherical 
basis to contain four major shells beyond those in the tri- 
axial basis. The convergence of the truncation has been 
checked by comparing ^Z n y(2j + l)v{nlj), Y, n lj Jx 

(nlj) 

and Y,ni 3 -n'i'j' Qw(nlj;n'l'j') with (N), (J x ) and (Q 20 ) 
respectively. The differences are typically of the order of 
0.01%. 

Energetics. In Fig. 1 the SM, CHFB and experimen- 
tal gamma ray energies Ej(J) = E(J) — E(J — 2) are 
plotted as a function of the angular momentum J. The 
SM results nicely reproduce the experiment including the 
backbending seen at J = 10. On the other hand, the 
mean field values of follow the same trend as the ex- 
perimental and SM ones but they are shifted downwards 
by w 0.5 MeV. This means that the mean field dynamic 
moment of inertia (J'^(J) = 4/A_E 7 ) is similar to the 
SM and experimental ones although the static moment 
of inertia {J^ x \j) = (2 J — 1)/E 1 ) is on the average a 
factor 1.5 bigger. (The origin of this discrepancy will be 
explained later.) 

Quadrupole properties. The striking similarity be- 
tween the SM and CHFB results up to the backbend can 
be gathered from the lower part of fig. 2, in which the 
intrinsic quadrupole moment is plotted along the yrast 
band. The SM values are extracted from the BE2 val- 
ues, assuming K = 0. The existence of an intrinsic state 
common to the members of the band can be guessed di- 
rectly by calculating the contribution of a given configu- 
ration to each SM wave function, {i.e., by summing the 



square of the amplitudes of all basic states having the 
same number of particles in each subshell). These con- 
tributions are practically identical in all the eigenstates 
up to J = 10. At higher spins rapid changes occur, and 
the configuration in which all the particles are in the fa / 2 
orbit becomes increasingly dominant. It is clear that the 
intrinsic state is becoming J-dependent at the backbend- 
ing region, and the discrepancies in fig. 2 beyond J = 10, 
suggest that it is no longer possible to extract an intrinsic 
Qo from the SM results assuming a K = band. In the 
upper part of the figure an alternative is proposed, by 
comparing the B(E2) values, obtained directly in the SM 
case with those derived from CHFB by applying the gen- 
eralization of the rotational model prescription to small 
triaxialities (see || ) . The agreement is again nearly per- 
fect up to J = 10 but then deteriorates, although not as 
much as in the lower figure. 

In assessing the significance of these results we should 
keep in mind that they are in both cases (rotational) 
model dependent. They indicate that the model is as 
good as exact up to the backbend, and then breaks down 
- at least in the standard implementation proposed here. 
They certainly do not indicate that the SM and CHFB 
descriptions are becoming different. On the contrary we 
shall find evidence of their closeness. 

Orbital occupancies. In figure 3 are have plotted 
the fractional occupancies of the spherical orbits in the 
CHFB solution (eq. (g )) (upper part) and in the SM 
one (lower part). In all cases they are quite constant up 
to the backbending where the fa / 2 orbit becomes rapidly 
the only relevant one. 

However: the fa/ 2 occupancy is always the largest by 
far , and in the CHFB case the contribution j x (f 7/2) to 
{J x ) in eq. (|3j) is always greater than 99 %. It means 
that the fa/2 orbit plays a major part in the two yrast 
regimes: below backbend as the major contributor to 
the deformed wavef unctions; and above through the fy 2 
configuration that becomes increasingly dominant. This 
picture is consistent with the usual idea that the back- 
bend is associated with alignment of fa/2 particles, which 
are also massively present in the collective regime at low 
spin. 

Magnetic properties. In Figure 4 we present the 
CHFB and SM results for the gyromagnetic factor g. In 
both cases and up to the backbending zone they are close 
to the rotational limit gn = Z/A = 0.50. For a pure fa/2 
configuration the value of g is also constant and equal to 
0.55 explaining the slight increase in g as we enter the 
backbending region where these configurations become 
dominant. 

Pairing properties. From all we have said, it follows 
that the SM and CHFB results are basically the same, 
except for a difference in the static moment of inertia . 
Its origin can be understood by redoing the SM calcu- 
lations reducing the JT = 01 two-body matrix elements 
involving orbits r and t according to 
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W% tt + 0.165V(> + l/2)( Jt + l/2), (5) 



02o/Q 2 oP/(HF) = 1.70(J - 0) • • ■ 1.63(J = 14). 



which amounts to subtracting a standard pairing term 
(j r is the angular momentum of orbit r). The resulting 
E 1 pattern for an exact calculation with the modified 
interaction is shown as SM(E) in fig. 5. To gain further 
insight we have also calculated in SM(P) the energies by 
taking expectation values of the modified interaction (|5|) 
using the SM wavefunctions obtained with KB3. (The 
coefficient 0.165 was chosen - somewhat arbitrarily - to 
make the first point coincide for CHFB and SM(P)). The 
conclusion is that: 

Although the energetics of the yrast band are strongly 
affected by the pairing modifications, the other proper- 
ties are not, since the wavefunctions change little. (The 
overlaps < SM(E), J\ SM, J > exceed 0.97 in all cases). 

The large static moments of inertia obtained in the 
CHFB calculations should be attributed to an inadequate 
treatment of pairing effects in a weak correlation regime: 
Exploratory tests using the Lipkin-Nogami approach on 
top of the CHFB scheme suggest that it is not the Gogny 
force that is responsible for the discrepancies but the lim- 
itations of the mean field treatment. 

48 Cr as axial rotor It has been recently argued JlJ] 
that the building blocks of wavefunctions describing good 
rotors are constructed by allowing particles to move in 
spaces defined by Aj — 2 sequences of major shell orbits, 
starting on the one with the largest j. For these blocks, 
an approximate form of SU3 symmetry is valid (quasi- 
SU3). One of the predictions of this model is that 48 Cr is 
an axially symmetric rotor, contrary to what happens to 
its counterpart in the sd shell, 24 Mg, that obeys Elliott's 
SU3 and is triaxial. Experimentally no second 2 + state 
is found in 48 Cr at low excitation energy, while in 24 Mg 
the second 2 + is degenerate with the yrast 4 + . In figure 
6 we present the values of the deformation parameters (3 
and 7 coming from the CHFB calculation. 

At first, [3 stays constant at (3 « 0.3, while 7 w 
which means that 48 Cr behaves indeed as an axial rotor 
up to the backbend. Above it, as f3 decreases fast and the 
system moves to a spherical regime making it difficult to 
interpret in a simple way the 7 behaviour. 

Effective charges. Finally, we can separate from the 
total quadrupole moment Q20 in CHFB, the valence con- 
tribution Q2oP/(HO) by summing q2o{n, I, j; n' , V , j') in 
eq. (||) over the 0/ and lp orbits, i.e., by identifying the 
valence orbits with harmonic oscillator ones. The ratio 

Q 20 /Q 2 oP/(HO) = 1.99(J - 0) • ■ ■ 1.83(J = 14) 

is quite consistent with the isoscalar effective charge used 
in the SM calculations q v + q^ = 2. 

Alternatively, we can define Q2op/(HF) by summing 
over all the values of I, j and j' corresponding to the 
pf shell, which amounts to use spherical HF orbits. This 
choice naturally reduces the effective charges but they 
remain quite constant since 



Conclusions. The quantitative equivalence of the SM 
and CHFB descriptions has two direct and welcome con- 
sequences: 

• It suggests that the Gogny force must be reasonably 
close to the realistic ones, consistent with NN data 
and known to yield high quality spectroscopy once 
their bad monopole properties are corrected. 

• It confirms the validity of the SM choice of a model 
space restricted to orbits in the vicinity of the Fermi 
level. 

Clearly there is much to be gained by combining the sim- 
plicity and rigour of CHFB with the SM precision and 
generality. 
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FIG. 1. Yrast energies £ 7 = E(J) - E(J - 2). 
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FIG. 2. Comparing B(E2) and Q trends. 
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FIG. 4. Gyromagnetic ratios 
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FIG. 5. Influence of the pairing strength on the moment of 
inertia. See text. 
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FIG. 3. Orbital occupancies. 



FIG. 6. The CHFB deformation parameters. 
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